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Buoyancy produced by autocatalytic reaction fronts can produce fluid flows that advect the front 
position, giving rise to interesting feedback between chemical and hydrodynamic effects. In a large 
diameter, extended cylinder that is relatively free of boundary constraints, localized initiation of an 
iodate-arsenous acid (IAA) reaction front on the bottom boundary generates a rising autocatalytic 
plume. Such plumes have several differences from their non-reactive counterparts. Using numerical 
simulation, we have found that if reaction is initiated using a spherical ball of product solution well 
above the bottom boundary, the subsequent flow can evolve much like an autocatalytic plume: the 
ball develops a reacting head and tail that is akin to the head and conduit of an autocatalytic plume, 
except that the tail is disconnected from the boundary. In the limit of large initial autocatalytic 
balls, however, growth of a reacting tail is suppressed and the resemblance to plumes disappears. 
Conversely, very small balls of product solution fail to initiate sustained fronts and eventually 
disappear. 



PACS numbers: 47.20.Bp, 47.70.Fw, 47.15.-x 

When feedback between autocatalytic chemi- 
cal reaction and fluid flow occurs, the result- 
ing "chemo-hydrodynamics" can lead to novel 
flow phenomena and interesting instabilities. The 
buoyancy produced by an autocatalytic reaction 
can drive fluid flows that deform the reaction 
front in an otherwise quiescent fluid. An extreme 
example of this deformation is the formation of 
rising plumes with complex morphologies and dy- 
namics. In this paper, we briefly review work 
on fluid flows driven by the iodate-arsenous acid 
(IAA) reaction, an archetypal cubic autocatalytic 
reaction. We consider various geometries with an 
emphasis on the rich dynamics of autocatalytic 
plumes. We then present numerical simulations 
of the effect of buoyancy driven flow on initially 
spherical "autocatalytic balls". Evolution of the 
flow structure resulting from an autocatalytic ball 
exhibits different morphological regimes that de- 
pend on the initial size of the ball. Extremely 
small initial balls die away and fail to produce 
sustained reaction fronts or fluid flow. Over a 
certain range of initial radius, the autocatalytic 
ball grows to resemble the structure of an auto- 
catalytic plume, with a well-defined head followed 
by an elongating tail. Beyond this regime, for a 
larger initial ball radius, a prolonged tail does not 
form, and the majority of reaction product re- 
mains confined within an ascending vortex ring. 



I. INTRODUCTION 

The propagating front in the IAA reaction consumes 
solution that is in a state where essentially no reaction 
has occurred, and leaves in its wake a nearly fully reacted 
solution that is less dense than the reactant. Heat is pro- 



duced during the reaction, as it is slightly exothermic. 
In addition, there is also an isothermal density change 
due to the difference in partial molal volumes of the re- 
actant and product solutions. In the presence of gravity, 
a front configuration where less dense product solution 
is beneath denser reactant solution is hydrodynamically 
unstable to buoyancy driven convection. Convection typ- 
ically deforms the front and causes it to propagate more 
quickly than it would in the absence of flow. 

The buoyancy-enhanced propagation of the IAA front 
has been carefully studied in various geometries. In thin 
vertical tubes, the relative strengths of the buoyancy and 
viscous forces due to the confining walls of the tube con- 
trol the stability of an ascending front. The front shape 
can be flat, asymmetric, or axisymmetric with respect to 
the tube axis [TH3]- When buoyancy is suppressed by a 
small tube radius, the front remains flat. At tube radii 
just above the onset of convection, the front becomes 
asymmetric. Further increasing the tube radius causes 
the front to develop an axisymmetric deformed shape, as 
shown in Fig [T] where the corresponding fluid flow is a 
vortex ring spanning the tube that moves with the front. 
We used this axisymmetric state as a test case for our 
numerical scheme, discussed in Section IV below. 

Regardless of their morphology, all ascending fronts in 
tubes have the property that their shape remains con- 
stant as they propagate. Since the shape is constant, 
front propagation is essentially a one-dimensional mo- 
tion, where the velocity of the front only depends on the 
relative buoyancy, the tube radius, and the reactant con- 
centrations E] . When a front is propagated against 
an imposed Poiseulle flow, its speed and shape also de- 
pend on the direction and magnitude of the imposed 
flow [5]. Surprisingly, in the wide gap limit for Poiseulle 
flow, there is almost no effect on the traveling velocity 
when reacting fronts travel in the opposite direction of 
the imposed flow [SJ [7] . 
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FIG. 1: (color online) a) The axisymmetric dimensionless 
concentration field is reflected across the axis to show the 
shape of the reaction front for a tube with radius 20£. b) The 
dimensionless vorticity field for the same front. Negative vor- 
ticity corresponds to clockwise rotation. A buoyancy- driven 
axisymmetric vortex ring is found in the vicinity of the front, 
with the fluid rising on the axis. The simulation used to gen- 
erate this image is discussed in Section IV. I is a length scale 
which is on the order of the front thickness. 



In the quasi two-dimensional limit of thin vertical slots 
or horizontal layers, a convecting reaction front takes 
on a two-dimensional pattern. As in the case of thin 
tubes, front motion and morphology in a thin horizon- 
tal layer depends on the chemical concentrations of the 
reactants [H HI IS]- Simulations show that for a com- 
pletely confined horizontal layer of solution, the propa- 
gation speed and deformation of the front depends on 
layer thickness pj]. In simulations of solutions with a 
free surface, surface tension gradients across the front 
can induce capillary (Marangoni) flows [IT] . 

Rising reaction fronts that span the width of a thin 
vertical slot (a Hele-Shaw cell) are unstable and will de- 
velop a wave-like pattern of fingers [T2TU4) . The fingers 
emerge from an initially flat front through competition 
between buoyancy-driven instability, which acts to ex- 
tend the fingers vertically, and diffusion, which smoothes 
the pattern [T5HTB] . Scaling of the fingers depends on 
the dimensions of the slot and reaction parameters [Hj . 
Amplification of the finger pattern can be achieved by 
matching the spacing of the the fastest growing mode 
with periodic variations in slot thickness [20] . Interest- 
ingly, descending fronts in a vertical slot can also have 
buoyancy-driven instabilities [3TJ[55], even though the de- 
scending buoyant front is initially stably stratified. The 
pattern results from the difference between molecular and 
thermal diffusivities. 

In the quasi one- and two-dimensional scenarios, mo- 
tion of the reaction front is severely limited by the bound- 
aries of the reaction geometry. We now focus on reaction- 
driven flows that develop in a three-dimensional geom- 
etry where boundary constraint is minimal. We first 
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FIG. 2: An autocatalytic plume produced by the IAA re- 
action. Starch indicator was used in solution to reveal the 
plume of product solution. 



discuss the behavior of autocatalytic plumes, and then 
compare them to simulations of a closely related system: 
autocatalytic balls. 



II. AUTOCATALYTIC PLUMES 

From the flow of air above a burning match to the ris- 
ing molten rock that forms volcanic island chains, plumes 
at vastly different scales are abundant in Nature. Under 
laminar flow conditions, plumes consist of a well-defined 
plume head and a trailing conduit, as shown in Fig. [2] 
Plume formation typically results from continuous forc- 
ing provided by a localized source of buoyancy. Usually, 
buoyancy is the result of a density difference caused by 
the thermal expansion of the fluid. Similarly, a source of 
fluid with a different chemical concentration will produce 
a plume due to the isothermal density difference of the 
fluid. Autocatalytic chemical reactions, such as the IAA 
reaction, produce both thermal and concentration den- 
sity differences, and therefore can induce complex plume 
formation. 

The upwelling flow from a laminar plume conduit pro- 
vides a continuous source of buoyancy to its head. It is 
usually the case that a plume head will remain attached 
to its conduit [23l El] , which sustains its ascent and feeds 
a vortex ring flow within it. If a vortex ring detaches from 
its source and becomes "free" , it is said to pinch-off. Free 
vortex rings are usually studied using a piston-cylinder 
arrangement [251 ES] , where a single piston stroke is used 
to create the vortical flow. However, free vortex rings can 
also be created by the pinch-off of autocatalytic plume 



3 



heads from their conduit [57J H5J. This process represents 
the most dynamically interesting stage of the evolution 
of an autocatalytic plume. These stages were described 
in detail in Refs. [57J [55J , and are summarized below. 

Autocatalytic plumes can be formed by initiating the 
IAA reaction at the bottom of a vertically oriented cap- 
illary tube and letting the reaction front escape into the 
bottom of a much larger cylindrical tank containing a 
quiescent volume of reactant solution. Plume formation 
is initiated once the reaction front emerges from the sub- 
merged capillary tube. The earliest stage of plume evolu- 
tion is by far the slowest, as the creeping flow of the tiny 
plume emerging from the capillary tube develops a small, 
roughly spherical head. The next stage involves the de- 
velopment of vortical motion in the head, where the up- 
welling product solution wraps itself toward the centre of 
the overturning vortex ring. At this point, the fate of the 
plume head — whether or not it pinches off, and if so, 
how many times pinch-off is repeated — depends on the 
fluid properties of the reactant solution. In cases where 
pinch-off occurs, the product solution flowing upwards in 
the conduit left behind forms a second generation plume 
head at the site where the conduit pinched. 

The evolution of autocatalytic plumes has been stud- 
ied in fluids with various properties by altering the con- 
centration of glycerol in the reactant solution [55]. In 
general, pinch-off does not occur for solutions with low 
concentrations of glycerol. These solutions have the least 
viscosity and the fastest diffusion constants. At the other 
limit, pinch-off was seldom observed for the highest glyc- 
erol concentrations. However, pinch-off does occur reli- 
ably in an intermediate viscosity range for glycerol con- 
centrations of 30% and 40% by volume. The highest 
number of pinch-offs observed for a given run was 5, for 
a 40% glycerol solution. The number of pinch-offs ob- 
served during a given experimental run increased with 
the height of reactant fluid for a cylindrical tank with 
the same diameter, reflecting the geometric limitations 
of even very large experimental volumes. 

Although the external appearances of autocatalytic 
and non-reacting laminar plumes are similar, there are 
many important differences between them. A comparison 
of the steady-states reached by plume conduits was pre- 
sented in Ref. [29J . Autocatalytic plume conduits have 
reaction occurring along the entire interface between re- 
acted and unreacted fluids, which causes the buoyancy 
flux to be distributed along the entire conduit length. 
This results in a velocity profile where the maximum ve- 
locity is near the interface between reacted and unreacted 
fluids. This contrasts the conduit velocity profile for a 
non-reactive plume, where the a maximum velocity is lo- 
cated along its axis of symmetry. The other major differ- 
ence between the two types of conduits are their steady- 
state morphologies. Due to the ongoing reaction along 
the front that forms the conduit, autocatalytic plumes 
have conical conduits that widen with distance from the 
base of the plume, whereas non-reacting laminar plume 
conduits are cylindrical. 



During the transient phase of plume growth, a non- 
reacting plume head is supplied with buoyancy by the 
upward flow in the attached conduit. This flow must 
rise at a greater velocity than the upward speed of the 
head. The plume head of a continuously driven, non- 
reacting plume is therefore unable to detach from its con- 
duit. Moreover, after a short initialization period, the up- 
ward velocity of a non-reacting plume head is constant, 
and its dimensions maintain the same scaling relation- 
ship as it grows [24]. On the other hand, autocatalytic 
plume heads accelerate, pinch-off, and do not maintain 
self-similarity [55]. The events leading to autocatalytic 
plume pinch-off were investigated in detail using simula- 
tion, and it was shown that a minimum velocity in the 
conduit beneath the head signals the onset of the pinch- 
ing off process long before morphological indicators, such 
as narrowing of the conduit |28| . Despite development of 
this velocity minimum, however, the head continues to 
accelerate upwards, demonstrating the lack of influence 
that the conduit has on the head at this stage of auto- 
catalytic plume development. The "conduit" therefore 
no longer acts as a conduit: it does not transport buoy- 
ant fluid to the head. Instead, the fluid trailing the head 
is more like a "tail" being dragged along in its wake. In 
Section V.B below, we numerically explore the growth of 
such tails on "autocatalytic balls" — buoyant, ascend- 
ing spherical regions of fluid bounded by autocatalytic 
reaction fronts. 



III. AUTOCATALYTIC BALLS 

The simplest geometry for a three-dimensional reac- 
tion front is a sphere. In the absence of buoyancy-driven 
flow, spherical autocatalytic reaction fronts have previ- 
ously been studied both theoretically and numerically. 
These studies have focused on the existence and stabil- 
ity of spherically symmetric solutions of the reaction- 
diffusion equations [501 |3"T] . and on the effect of auto- 
catalyst decay , 32J . Threshold conditions needed to form 
spherical autocatalytic reaction fronts have also been cal- 
culated for both quadratic and cubic autocatalytic reac- 
tion fronts [3"3"] . 

Autocatalytic reaction fronts under gravity have some 
similarities to flame fronts [MJ [35] . Both systems involve 
positive feedback mechanisms that sustain chemical reac- 
tions occurring at their respective fronts. In the presence 
of an oxidizing agent, heat produced by exothermic chem- 
ical reactions in a combustible system sustains the reac- 
tion process. This process is analogous to the production 
of catalyst by an autocatalytic reaction front, which pro- 
duces the catalyst needed to maintain propagation of its 
reaction front. As a consequence of this analogy, spher- 
ical autocatalytic reaction fronts are sometimes referred 
to as "isothermal flame balls" [301 [35] . They have also 
been called "spherical reaction balls" [3T|. In this paper 
we adopt the name autocatalytic balls. 

Flame balls are steady, radially symmetric solutions 
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of the reaction-diffusion-conduction equations for pre- 
mixed laminar flames. Originally proposed by Zel- 
dovich in 1944 36 , flame balls have since been observed 
when a spherical flame is formed under free fall condi- 
tions [23 [35J. The combustive process which sustains 
flame balls is highly exothermic, and therefore free fall 
conditions are necessary to prevent vigorous buoyancy 
driven flow. During flame ball experiments in micrograv- 
ity, small local accelerations create unwanted gravity-like 
random accelerations called "g-jitter". In the presence 
of these small acceleration fluctuations, flame balls are 
sometimes observed to deform into "flame strings" [3T1 - 
131)] . Flame strings are simply flame balls that have been 
stretched into a long, cylindrical shape. 

In addition to being caused by unwanted accelerations 
in microgravity experiments, the deformation of flame 
balls under gravity is particularly relevant to explosion 
scenarios for Type la supernovae 40J . The initial stages 
of the march toward detonation of a Type la supernova 
involve the formation of a "flame bubble" near the centre 
of a compact white dwarf star. This reacting bubble 
is driven by its buoyancy toward the nearest surface of 
the star in the early stages of the supernova explosion, 
breaking its spherical symmetry. 

IV. AUTOCATALYTIC BALL SIMULATION 

To simulate the behavior of autocatalytic balls, we em- 
ployed a dimensionless model previously used to study 
autocatalytic plumes [28]. The model uses the semi- 
implicit method for pressure linked equations (SIMPLE) 
algorithm [1TJ [12] to calculate the relevant physical fields. 
The formulation of the model is detailed in this section. 

The dimensionless parameters required for simulating 
autocatalytic balls come directly from the scaling of the 
basic dynamical equations. The thermal, chemical, and 
mechanical properties of the fluid are characterized by 
the Schmidt number 

Sc = v/D , (1) 

and the Lewis number 



These values of Sc and Le indicate that the rates of dif- 
fusion of heat and momentum are much faster than that 
of concentration, implying that the concentration profile 
at the reaction front has by far the steepest gradient of 
the relevant fields. 

The Boussinesq approximation was used to quantify 
the density change from the reaction as 

p = p [l-aAT-pAc], (5) 

where AT = T — Tq and Ac = c — Co are the temperature 
and concentration differences across the reaction front, 
and po, Tq, and Co are the initial conditions for the den- 
sity, temperature and concentration, respectively. Here, 
a is the thermal expansion coefficient and /3 is the com- 
positional expansion coefficient. Since the density change 
comes from a linear combination of temperature and con- 
centration effects, these contributions can be separated 
into two Rayleigh-like numbers [2T] [22] 

Ra T = ^-AT and Ra c = ^-Ac (6) 
vD vD w 

where I = \J~Dt is the molecular diffusion length scale of 
the front, and r is a time scale that depends on the front 
velocity Vf, such that Vf = y/D/r. The simulations in 
this paper used 

Ra T = 0.10 and Ra c = 0.17 , (7) 

which were calculated using the parameter values for an 
IAA solution containing 40% glycerol given in Ref. [25] . 

Autocatalytic balls were assumed to be axisymmctric 
about the direction of gravity, allowing the simulation to 
be carried out in the cylindrical coordinates (r, z), where 
gravity is oriented in the negative z direction. In this co- 
ordinate system, the set of dimensionless equations used 
to model autocatalytic ball evolution were 

V • u = 0, (8) 
^-* = -Vp + \7 2 u+{Ra T T + Ra c c)e z , (9) 



Le = k/D , (2) 

where v is the kinematic viscosity, k is the thermal diffu- 
sivity, and D is the molecular diffusivity of the autocat- 
alyst, which for the IAA reaction is the iodide ion. The 
ratio of Sc and Le is the Prandtl number, 

Pr = v/k , (3) 

which is a useful quantity for characterizing purely ther- 
mal plumes [13] . 

For a 40% glycerol IAA solution, using k = 1.2 x 
10~ 3 cm 2 /s and values from Ref. [25], we find 

Sc = 9000, Le = 280, and Pr = 32. (4) 



— = V 2 c + ^(c), (10) 

and 

DT 

— =LcV 2 T + F( C ). (11) 

Here, e 2 is a unit vector pointing in the z direction and 
D/Dt is the material derivative. The cubic autocatalytic 
reaction term is 

F( C ) = c 2 (l-c) , (12) 

where the dimensionless concentration c is normalized 
by the maximum possible concentration of the product 
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FIG. 3: (color online) a) For an autocatalytic ball with ro = 
21, the concentration profiles at different times, in units of r, 
and b) the change in the maximum value of the concentration 
c m ax with time, c) The concentration profiles at different 
times for an autocatalytic ball with ro = 3£. Red curves 
indicate concentration profiles when c max is decreasing, blue 
curves are profiles when c m ax is increasing, d) The change in 
Cmax over time for the ro = 3£ ball. 



solution. This normalization scheme gives a maximum 
dimensionless concentration of c = 1 and a minimum of 
c = 0. The dimensionless temperature T is normalized 
in a similar way as the concentration, and is given by 
scaling the physical AT by ATjj = —AHAc/poC p , which 
bounds it by < T < 1. Here, AH is the heat of reaction, 
and C p is the specific heat at constant pressure. 

The cylindrical boundaries of the system are given by 
r = rh, z = and z = Zb- Conditions at the bound- 
aries are defined by imposing no flow of concentration or 
temperature across the boundary, and no-slip, giving 



u = 



u 



dT dc n 

— = — = at r = r b , 

or or 



(13) 



dT 

dz 



dc 
dz 



= at z = and z = 2%. (14) 



The rather small length £, which is of order the thick- 
ness of the reaction front, is used to scale all of the other 
lengths in the model. This means that rather large di- 
mensionless values of r^ and 2& are required to model an 
extended volume of fluid. 



V. SIMULATION RESULTS 
A. Front death 

Small localized regions of autocatalyst are not capable 
of sustaining three-dimensional reaction fronts for an au- 



tocatalytic reaction with cubic kinetics [33]. This "front 
death" phenomenon is caused by diffusion, and is un- 
like the threshold behavior of excitable systems, such as 
the action potential of neurons |H] or the propagating 
oxidation wave in the Belousov-Zhabotinsky (BZ) reac- 
tion [45]. Neurons and the BZ reaction are excitable in 
the sense that they remain in a stable state until a thresh- 
old is reached, sending them into an unstable state. For 
a neuron, excitation occurs once a certain voltage, or ac- 
tion potential, is reached, resulting in the transmission of 
an electric signal. The excited state of the BZ reaction 
is reached when enough catalyst is present to transmit 
a propagating oxidation pulse through the solution. A 
BZ reaction pulse does not consume all of the reactants 
as it moves through the reactant medium, allowing for 
multiple excitations. A cubic autocatalytic reaction like 
the 1AA reaction, in contrast, always has locally unsta- 
ble kinetics [33] . Therefore the presence of any amount 
of catalyst will render the solution unstable. However, 
despite the localized instability created by the presence 
of autocatalyst, three-dimensional diffusion can prevent 
the initiation of a reaction front provided that a critical 
amount of catalyst has not been exceeded [IB]. In this 
section, the consequences of front death on autocatalytic 
balls is explored using our numerical model. 

To initiate autocatalytic balls, spheres located on the 
axis of symmetry of the coordinate system with c = 1 
and various initial radii r were used as initial condi- 
tions. Convective effects were suppressed in the simula- 
tions exploring front death, so that it could be studied as 
a pure reaction-diffusion phenomenon. Even when buoy- 
ancy is introduced into the calculation, the buoyant force 
is negligible on small spheres in the radius range where 
the front dies, and very little convective transport occurs 
during the time scale it takes for the front to die. All 
simulations investigating front death were performed on 
a spatial domain of 2% = r& = 120£, with the centre of the 
initial ball located on the axis of symmetry at zq = 60£. 

The radial distribution of product solution at different 
times for an autocatalytic ball with ro = 2£ is shown in 
Fig. [3^,. The time evolution shows that the product con- 
centration simply diffuses away from its initial position 
without forming a reaction front. The maximum value 
of the concentration, c maxi is shown for the same ball in 
Fig. [3]d. The decrease in c max towards c = indicates 
death of the reaction front. An autocatalytic ball with 
an initial radius of ro = 3£ has a different fate. In this 
case, the initial front begins to diffuse away as it did for 
r = 11. After a period of time, however, a sustained 
reaction front is initiated and the autocatalytic ball re- 
covers from the initial concentration drop. This is seen 
in the radial distribution of product concentration for 
various times shown in Fig. [3]:, and the evolution of the 
maximum concentration shown in Fig. [3ji. At times less 
than ~ lOr the concentration diffuses and insufficient 
new product is generated. After that time, however, the 
front recovers, leading to the stable outward propagating 
front structure shown in Fig. [3fc at 40r and 50r. 



6 



2000 




FIG. 4: The a) concentration, b) temperature, and c) vortic- 
ity fields at time 20r for an ascending autocatalytic ball with 
ro = 151. The spatial domain of the simulation was r b = 500£ 
and Zb = 2000£, however only data within the coordinates 
r < 1001 and 200£ < z < 6001 is shown. Concentration and 
temperature data are shown reflected across the z axis, vor- 
ticity data, where negative values indicate clockwise motion, 
is not reflected. 



B. Heads and tails 

When buoyancy is suppressed and front death does not 
occur, autocatalytic balls maintain a spherical shape as 
they expand. We now turn our attention to the evolution 
of autocatalytic balls under gravity. This involves an 
interesting dynamic interplay between reaction, diffusion, 
and buoyancy-driven flow. For these simulations, only 
autocatalytic balls that do not die are considered. 

An example of the concentration, temperature, and 
vorticity fields from a simulation of an autocatalytic ball 
with ro = lhl are shown in Fig. |4j The large ratio of the 
diffusivities k and D (i.e. the large Lewis number Le) is 
evident in the difference between the concentration and 
temperature fields, the temperature field being the more 
diffuse of the two. The structure of the vorticity field in 
Fig. |4j; is typical of the ring vortex that develops in all 
ascending autocatalytic balls. 

The location of the top edge of the concentration jump 
for autocatalytic balls with different initial radii ro is 
shown in Fig. [5| The centres of the balls in these sim- 
ulations were initially located at zq — 100^ in a spatial 
domain with ri, = hOOl and z b = 2000£. The simulation 
time-step was At = 0.005t. These are the same param- 
eters used in Ref. |28j to simulate autocatalytic plumes, 
with the exception of the initial ball location z , which 
was farther from the lower boundary of the computa- 
tional domain than the one used for plumes. The auto- 
catalytic ball with the largest initial radius ro reaches the 
top of the computational domain faster than the others. 
The smaller the ro of the ball, the longer it takes to reach 
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FIG. 5: The ascent of autocatalytic balls with different initial 
radii over time. The curves correspond to autocatalytic balls 
with initial radii r of a) 27.5^, b) 22. 51, c) 17. 5£, d) 15£, e) 
10£, f) 7.51, and g) 51. 



the top. All of the balls reach a terminal velocity that 
is approximately the same, increasing only slightly with 
ro. The range of terminal velocities was 64.9 to 68.9 in 
units of l/r, which corresponds to a physical velocity dif- 
ference of only 5 x 1CP 3 cm/s between the smallest and 
largest ro. The terminal velocities are dictated by the 
final amount of product solution in the ascending flow 
structure. As the small range of terminal velocities sug- 
gests, the amount of reacted fluid in the vortex ring once 
it reaches the top boundary is similar for all values of 
ro. The integral of the concentration field for the vortex 
ring indicates that the amount of product reaching the 
top for the ro = hi and ro = 27. hi balls was different by 
only ~ 8%, with the larger ball delivering a slightly larger 
amount of product. Even though a larger ro balls starts 
with a larger amount of product, the slow initial ascent 
of the smaller ones gives them more time to expand and 
accelerate before reaching terminal velocity of their way 
to the top. 

The evolution of the morphology of a small autocat- 
alytic ball with ro = hi is shown in Fig. [6] In the initial 
stages of its evolution, it remains roughly spherical. How- 
ever, as it ascends it leaves behind a tail: a thin trail of 
product solution extending into the wake of the rising au- 
tocatalytic ball. As the tail elongates, the upper portion 
of the ball develops into a vortex ring that resembles a 
plume head. Just as an autocatalytic plume head pinches 
off from its conduit, the autocatalytic ball head pinches 
off from its tail. As it does so, the top of the tail grows 
a new, second-generation head, and the first-generation 
head rises upwards as an essentially free vortex ring. 

The similarity between the morphological evolution of 
autocatalytic plumes and balls disappears for balls with 
larger initial radii. A ball with r = 27 .hi, which fails to 
develop a reacting tail, is shown in Fig. [7} As it rises, the 




FIG. 6: The evolution of an autocatalytic ball with ro = 5£, 
showing the development of the head and a tail of reacting 
solution. The spatial domain of the simulation is — 5Q0£ 
and Zb = 200(M, however the portion of the field displayed 
field shown has a radius of 1001 From the left-most frame to 
the right-most frame, the elapsed time after initiation is lOr 
to 70r in increments of lOr. 



FIG. 7: The evolution of an autocatalytic ball with ro = 
27.5£, showing that only a very thin filament of reacted so- 
lution trails the head, and that the filament is not capable 
of producing a sustained reaction front. The spatial domain 
of the simulation is rb = 500^ and Zb = 20001 however the 
portion of the field displayed has a radius of 1001 From the 
left-most frame to the right-most frame, the elapsed time after 
initiation is 2.5r to 32. 5r in increments of 5r. 



ball leaves behind a filament of reacted solution. How- 
ever, unlike the tails formed in the wake of autocatalytic 
balls with r < 25£, the filament is so thin and dilute 
that it diffuses away before it can initiate a reacting front, 
much like how a very small autocatalytic ball undergoes 
front death. While the exact cutoff radius for reacting 
tail development, rot, was not determined, the simulation 
results show that it lies in the range 25^? < rot < 27. 51 

To help understand why tail production ceases once 
ro surpasses rot , it is useful to use a dimensional analysis 
argument for autocatalytic balls similar to one developed 
for astrophysical flame bubbles [10] ■ For highly viscous 
autocatalytic balls, two time scales can be compared to 
estimate the conditions when a spherical ball of prod- 
uct undergoes significant deformity from viscous effects. 
These time scales are the viscosity time scale 

To = ^, (15) 

v 

and the buoyancy time scale 




where g 1 = g(p u — Pr)/pu is the buoyancy of the flame 
ball, with g being the gravitational acceleration, and p r 
and p u being the densities of the reacted and unreacted 
fluid, respectively. The ratio of the two time scales 

n v -3/2 ,,„ A 

V = 77 r ° (17) 
* v v g 

shows which time scale dominates as a function of the 
initial radius ro. For the case where tj, <C t v , the ball 
travels upwards for a significant period of time before 
viscous effects to are able to alter its shape. The buoyant 
and viscous timescales are equivalent at a ball radius of 

r.-( F ) ■ (18) 

Using the parameters of the simulation, and a total 
density jump across the front of Ap = (p u — p r ) = 
6.4 x 1CT 4 g/cm 3 [47], we find r c - 331 This scale is 
in rough agreement with the observed crossover ro be- 
tween autocatalytic balls with reacting tails and those 
without. 
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tocatalytic balls was explored. Using reaction conditions 
for the iodate-arsenous acid reaction in 40% glycerol solu- 
tion, three different regimes of behavior for autocatalytic 
balls with different initial radii were found. Very small 
autocatalytic balls that undergo front death occupy the 
first regime. Buoyancy is negligible in this case. In the 
second regime, where 3£ < r < 25£, autocatalytic balls 
closely resemble autocatalytic plumes: they develop re- 
acting heads and tails. Pinch off can occur, producing 
secondary heads. In the third regime, for the largest ini- 
tial condition examined with r = 27.5^, only a small 
amount of product is dragged into the wake of the head, 
and the resulting filamentous structure does not produce 
a reacting tail. A dimensional analysis argument involv- 
ing the viscous and buoyancy time scales was used to 
help explain the disappearance of the tail. 



FIG. 8: The viscous time scale, t v , and the buoyancy time 
scale, Tb, as a function of autocatalytic ball radius ro. The 
shaded area is the region where the transition from balls with 
tails to those without is observed in the simulations. 



A plot of Tb and t v in the range of < ro < 401 
is shown in Fig. [8] For ro < r c , the timescales are such 
that t v < Tf,, and viscous forces have more time to act on, 
and therefore deform, the autocatalytic ball as it ascends. 
For r > r c , becomes shorter than t v and viscous 
forces have less influence on autocatalytic ball morphol- 
ogy before the buoyancy forces carry the ball upwards. 
In our simulations, viscous deformation of the autocat- 
alytic balls for r < 25£ was capable of dragging enough 
product out of the autocatalytic ball to form a tail with a 
sustained reaction front of its own. For ro > 27. 5£, how- 
ever, buoyancy forces dominate and viscous effects drag 
only a small amount of reacted solution into the wake of 
the rising ball. The resulting filament of product solu- 
tion diffuses away and is insufficient to initiate a reacting 
autocatalytic front. 



VI. CONCLUSIONS 

Using a numerical scheme previously employed to sim- 
ulate autocatalytic plumes, the evolution of buoyant au- 



The scope of our understanding of autocatalytic 
plumes and balls has thus far been limited to those 
formed by the iodate-arsenous acid reaction. It would 
be interesting to explore how these phenomena behave in 
other autocatalytic reactions, for instance in a reaction 
in which the thermal and concentration contributions to 
buoyancy have opposite signs. The resulting interplay 
between chemical reaction and fluid flow would almost 
certainly produce interesting and surprising dynamical 
phenomena. 
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